Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 5 February 2008 



(MN MfeK style file v2.2) 



Magnetically moderated outbursts of WZ Sagittae 



0. M. Matthews 1 *, R. Speith 2 , G. A. Wynn 3 and R. G. West 3 

1 Laboratory for Astrophysics, Paul Scherrer Institut, Wurenlingen and Villigen, CH-5232 Villigen PSI, Switzerland 

2 Institut fur Astronomie und Astrophysik, Universitdt Tubingen, Auf der Morgenstelle 10, D-72076 Tubingen, Germany 

3 Department of Physics & Astronomy, University of Leicester, University Road, Leicester, LEI 7RH 



5 February 2008 



O 

O 

(N 

> 
O 



> 

o 
o 

(N 



O 

S3 

Oh: 

6 

CO 



X 



ABSTRACT 

We argue that the quiescent value of the viscosity parameter of the accretion disc in 
WZ Sge may be a co id ~ 0.01, in agreement with estimates of a co id for other dwarf 
novae. Assuming the white dwarf in WZ Sge to be magnetic, we show that, in qui- 
escence, material close to the white dwarf can be propelled to larger radii, depleting 
the inner accretion disc. The propeller therefore has the effect of stabilizing the in- 
ner disc and allowing the outer disc to accumulate mass. The outbursts of WZ Sge 
are then regulated by the (magnetically determined) evolution of the surface density 
of the outer disc at a radius close to the tidal limit. Numerical models confirm that 
the recurrence time can be significantly extended in this way. The outbursts are ex- 
pected to be superoutbursts since the outer disc radius is forced to exceed the tidal 
(3:1 resonance) radius. The large, quiescent disc is expected to be massive, and to be 
able to supply the observed mass accretion rate during outburst. We predict that the 
long-term spin evolution of the white dwarf spin will involve a long cycle of spin up 
and spin down phases. 

Key words: accretion, accretion discs - binaries: close - stars: individual: WZ Sge - 
stars: magnetic fields. 



1 INTRODUCTION 

The SU UMa stars are a subclass of Dwarf Novae (DN) 
whose light curves exhibit superoutburst behaviour. Super- 
outbursts typically occur at intervals of several months, in- 
terspersed with normal outbursts every few weeks. In gen- 
eral, superoutbursts are around one magnitude brighter than 
normal outbursts and last a few weeks rather than 2-3 days. 
WZ Sagittae is an unusual SU UMa star because the recur- 
rence time is extremely long and no normal outbursts are 
observed between superoutbursts. Superoutbursts occurred 
in 1913, 1946, 1978 and 2001, with the last of these coming 
somewhat earlier than expected after the previously regular 
recurrence time t rcc ~ 33 years. 

Normal outbursts in DN can be explained by a thermal- 
viscous disc ins tability due to the partial ionisation of hy- 
drogen (see e.g. ICannizz o 1993, for a review). The standard 
disc instability model (DIM) assumes that the viscosity of 
the disc material increases by around one order of magnitude 
during outburst, enhancing the mass-fl ow rate through the 
disc. T he usual values adopted for the IShakura fc Sunvaevl 
1 19731 ) viscosity parameter are a co id ~ 0.01 in quiescence 
and ahot ~ 0.1 in outburst. In an extension to this the- 



ory, [Osaki (1995) proposed that superoutbursts are due to a 
thermal-tidal instability. During each normal outburst, only 
a small fraction of the total disc mass is deposited onto the 
white dwarf (WD), leading to an accumulation of disc mass 
and angular momentum. The outer disc radius therefore ex- 
pands until the tidal radius is reached. At this point, en- 
hanced tidal interaction with the secondary star is thought 
to increase dissipation in the disc and so raise the mass ac- 
cretion rate. A superoutburst may then be triggered. This 
model successfully explains the superoutburst behaviour of 
ordinary SU UMa stars using the same viscosity parameters 
as the standard disc instability model (DIM). 

ISmakl (1 19931 ) argued that the viscosity in WZ Sge 
must be far lower than the standard values (so that 
Qcoid < 5 x 10~ 5 ), for the following reasons: 

Smak 1 In the standard DIM, the inter-outburst re- 
currence time is governed by the viscous time-scale of the 
inner accretion disc, which is far shorter than the observed 
value of ~ 30 yr. Lowering a co id slows down the viscous 
evolution of the accretion disc, increasing the inter-outburst 
time. 
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Smak 2 Integrating the critical surface density using 
the standard DIM a co id gives a maximum disc mass that is 
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far lower than the observed estimate of the mass accreted 
during superoutburst, AAf acc ~ 10 24 g. Lowering the 
viscosity allows more mass to be accumulated in the outer 
regions of the accretion disc to fuel the outburst accretion 
rate. 



There is however no obvious reason why the quiescent 
accretion disc in WZ Sge should have a much lower viscosity 
than the other SU UMa stars. This has led a number of 
authors to suggest other reasons to account for the long 
recurrence time and high disc mass prior to superoutburst. 

In order to produce a long recurrence time without 
requiring a low a co id the inn er, most unstab le regio ns of 
the disc mus t be stabilized. iHameurv et al.l (|l997T j and 
IWarner et al.l (|l996l ) find that the disc can be stabilized if 
the inner regions are removed, either by evaporation into a 
coronal layer (see also lMever et al] ll999). or by the presence 
of a magnetosphere. This need not necessari ly lead to an in- 
crease d recurrence time however. Models bv lHameurv et all 
l| 19971 ) suggest a marginally stable disc that requires an 
episode of enhanced mass transfer to trigger a disc out- 
burst. The recurrence time in this model is governed by 
the mass transfer fluctuation cycle, which has no obvious 
physical connection with the regular 33 y ear outburst cycle 
observed in WZ Sge. IWarner et al.l l|l99rj ) also find that un- 
der certain conditions the disc will be marginally unstable 
and produce outside-in outbursts with the required recur- 
rence time (We explore the case of a marginally stable disc 
in Section [3.3p . However, neither of these models addresses 
the problem of how to accumulate enough disc mass during 
quiescence to explain the mass accreted in outburst (Smak 
2). Only ~ 10 21 g is available in the disc just prior to out- 
burst so that both of the above models require mass to be 
added to the disc during the outburst. The authors appeal to 
irradiation of the secondary star to increase th e transfer rate 
during outburst and supply the missing mass. IWarner et al.l 
(1996) only produce normal outbursts, not superout bursts 
in their model. The reason for this is that their disc radius 
J?disc ~ 1.1 x 10 10 cm is far smaller than the tidal radius, so 
enhanced tidal interaction with the secondary (required to 
produce a superoutburst) is not possible. 

In this paper we argue that the value of a co id in WZ 
Sge may be consistent with standard DIM values (a CO M ~ 
0.01). In se ction [5] we perform a similar calculation to 
Smakl (Il993l) adopting the system parameters suggested by 
Spruit fc Rutted (|l998h . For a fixed disc mass we find that 
Qcoid oc i?o U 9 t (where 72 out is the outer disc radius) , and show 
that a value of a co id « 0.01, and a moderate increase in -Rout 
results in a value for AM acc consistent with observation. We 
then obtain a similar value for a co id by estimating the cool 
viscous time-scale from the observed superoutburst p rofile. 

X-ray observations of WZ Sge bv lPattersonl (1998) have 
revealed a coherent 2 7.86 s oscillation in the ASCA 2-6 keV 
energy band (see also |Pattersonin~980j ) . The most convincing 
explanation for these data is that WZ Sge conta ins a mag- 
netic white dwarf. A paper bv lLasota et al.l (|l999l ) interprets 
WZ Sge as a DQ Her star and the 27.86 s oscillation as the 
spin period of the white dwarf (P S pm). The authors then 
suggest that WZ Sge is an ejector system: i.e. most of the 
material transferred from the second ary is ejected from the 
system, similar to the case of AE Aqr (|Wvnn et al.lll997h . In 



this model an accretion disc is recreated in outburst, which 
again is triggered by a mass transfer event. 

In section [3] we suggest that if the magnetic torque is 
not sufficient to eject mass completely from the system, but 
merely propels it further out into the Roche Lobe of the WD 
(we refer to the system as a weak magnetic propeller), then 
(Smak 1) and (Smak 2) can be explained with a standard 
DIM value for a C oid- The injection of angular momentum 
into the disc by the propeller radically alters the surface den- 
sity profile, compared with that of a conventional DIM disc. 
The inner disc in such a system is empty, and outbursts can- 
not therefore be triggered there in the normal way. The pro- 
peller also inhibits accretion, and mass accumulates in the 
outer disc, increasing the recurrence time significantly. The 
disc may also be forced to expand to the tidal radius (and 
a little beyond) which could allow further material to accu- 
mulate in the outer regions of the disc before triggering an 
outburst. This would increase the recurrence time and out- 
burst mass still further, solving (Smak 1) and (Smak 2). 
The outburst is initiated as in the standard DIM, and does 
not require any episode of enhanced mass transfer. Simula- 
tions are performed in one and two dimensions. The former 
are computationally cheap, so that we can present several 
outburst cycles, but cannot account for the enhanced mass 
storage due to tidally drive disc expansion. 

In Section [4] we discuss the observational evidence for 
a truncated accretion disc with reference to observed spin 
periods and spectral data. Both these results support are in 
agreement with a truncated disc. In section[5]we consider the 
spin evolution of the WD, and whether a spin cycle might 
exist. Finally, in section [6] the applications to other CVs are 
discussed. 

2 THE QUIESCENT VALUE OF a IN WZ SGE 

2.1 The quiescent disc mass as a limit on a co id 

Following ISmakl l| 19931 ) we estimate a value for a co id by con- 
sidering the disc mass immediately prior to outburst. An 
outburst is triggered when the surface density at some radius 
E (R) becomes greater than the maximum value (S cr it) al- 
lowed on the cool branch of the E —T relation (e.g. lCannizzol 
[1993J) at that radius. Therefore, the maximum possible disc 
mass prior to outburst is given by: 

/■flout 

M max ~ / 2nRT, clit (R)dR, (1) 

where R ln and i? ut are the inner and outer disc radii re- 
spec tively. Calculations of the vertical disc structure yield 
(e.g. iLudwig et a"i1ll994T ) 

V / r>\ C7n JJ- 0-37 ol'W -0.8 0.4 „-2 / n \ 

Ecrit (R) = 670 Mj R 10 a c>0 01 fj, g cm , (2) 

where Mi is the mass of the WD in solar masses, Rio = 
-R/10 10 cm, a c ,o.oi — a C oid/0.01, and /i is the mean molec- 
ular weight. Assuming that /i ~ 1 and setting Ri n = 0, 
equations JT} and @ can be combined to give 

M max ~ 1.36 x 10 23 Mf a37 a-S;oi i^ut.io • (3) 

Although WZ Sge has a long observational history, the val- 
ues of some of its fundamenta l parameters are still uncer- 
tain. For instance, [Smak (1993) deduces that Mi = 0.45 and 
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q = M2/M1 = 0.13, where as the spectroscopic observations 
of ISpruit fc Ruttenl (Il998|) (SR hereafter) le ad to values of 
Mi = 1.2 and q = 0.075. lLasota et alJ [| 19991 ) point out that 
a white dwarf mass of Mi = 0.45 is too small if it is assumed 
to be rotating with a period of 27.87 s (although Mi = 1.2 is 
possibly too high), so for the purposes of this paper we shall 
adopt the system masses suggested by SR. We also adopt 
an orbital period of 82 minutes. 

ISmakl l|l993l ) estimated the total mass accreted dur- 
ing outburst to be around AM acc ~ 1-2 x 10 24 g. Setting 
AM acc < M max , and using the system parameters Mi = 1.2, 
q = 0.075 and R out = 0.37a = 1.7 x 10 10 cm (SR), we find 
from equation © that a co id < 0.006. While this value is 
around a factor of 2 lower than the standard DIM value 
of a co id ~ 0.01, it is 2 orders of magnitude higher than the 
value obtained by Smak. However, setting M max ~ AM acc , 
we find « co id oc R ou t 3 ' 9 , i.e. a co id is very sensitive to disc 
radius. Therefore, any error in the observationally inferred 
value will produce a large discrepancy in a co id . The observa- 
tional estimate in SR (-Rout = 0.37a) was determined from 
models of the accretion stream disc impact region, which 
was found to be rather extended in the case of WZ Sge. A 
theoretical estimate of the tidal radius can be found by uti- 
lizing smoothed particle hydrodynamic calculations assum- 
ing the system parameters detailed above. This approach 
yields a tidal radius of -Rout — 0.5a (see section 1X4")) . Assum- 
ing this value for -R ou t and repeating the above calculation 
we find that a co id )$ 0.02, in good agreement with the stan- 
dard DIM valu e . We note that the analytic treatment of 
lOsaki fc Mever] (|2002h gives a similar R ou t when the same 
parameters are used. 



2.2 Extracting a co id from the outburst lightcurve 

A second, independent estimate of a C oid can be obtained us - 
ing the observed outburst profiles l|Patterson et al.l Il98ll ) . 
Superoutbursts have similar observational profil es to the 
outbursts of x-ray transients (|King fc Ri tter 1998), i.e. they 
have a rapid rise to outburst maximum, followed by an ex- 
ponential decay on the hot viscous time-scale until the disc 
is no longer hot enough to remain fully ionized and falls into 
quiescence. This produces a steep decline in brightness on a 
thermal time-scale followed by a slower decline as the disc 
readjusts on the co ol viscous t ime-sc ale 

Fig [2] (from iKuulkersI [2000) presents the visual 
lightcurves for the 1913, 1946 and 1978 outbursts. It is pos- 
sible to determine a co id by considering the visual luminosity 
as the disc returns to quiescence, i.e. ~ 30 — 90 days af- 
ter the outburst began. This is more difficult in the case of 
the 2001 outburst since the lightcurve is complicated by the 
presence of echo outbursts. The disc brightness decays expo- 
nentially on the cool viscous time-scale t v isc,c (this is clearly 
seen in the 1946 outburst, but is less obvious in 1978), such 
that the visual luminosity is given by: 



F (t) = F exp 



^visc,< 



(4) 



where t is time and the cold viscous time-scale is given by 

(5) 
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Figure 1. Schematic diagram of a superoutburst profile. 
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Figure 2. Plot of optical mag nitude dur i ng th e 1913, 1946 and 
1978 outbursts of WZ Sge, from lKuulkersl (|2000h . Reprinted with 
permission from Elsevier. 



The cold state viscosity is given by (Sh akura fc Sunvaevl 



a co i d c s H 



(6) 



Here c s is the adiabatic sound speed, and H is the disc 
scale height. During quiescence the disc temperature is T ~ 
5000 ± 1000 K. If we assume that fi ~ 1, the sound speed 

1 We consider the 



is then c s 



kT 



i.4 x 10 s cm s" 



part of the disc from which the majority of the emission 
should come during the cold viscous decay, namely R ~ 
2 x 10 ± 1 x 10 10 cm. Viscous dissipation is given by 

9 



D(R) 



(7) 



wher e 57 is the disc angular frequency (e.g. iFrank et al.l 
2002). The very inner part of the disc is evacuated, while 
the outer part has a low angular frequency so that our 
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Figure 3. Plots of dissipation, showing outburst behaviour, in 
one-dimensional simulations of the WZ Sge accretion disc. Frame 
(a) is without a stellar magnetic field and frame (b) frame has a 
field of B ~ 1 kG. 
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Figure 4. Plot of dissipation, showing a magnetically moderated 
outburst in one-dimensional simulations. 
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Figure 5. Plot showing evolution of the accretion disc mass as 
a function of time magnetic and non-magnetic cases. The mag- 
netic outbursts (where B ~ 1 kG) consume more than half of 
the accretion disc, about ten times as much as the non-magnetic 
outbursts. 



estimate for R seems sensible. We further assume that 
H/R ~ 0.10 ± 0.05 for a geometrically thin disc. From 
the 1946 outburst lightcurve it can be seen that, as the 
disc falls into quiescence, the visual flux drops by one mag- 
nitude in ~ 130 ± 40 days. From equation ([5]) we find 
Qcoid ~ 0.028 ± 0.022, in reasonable agreement with the es- 
timate in section [2711 



3 A MAGNETIC SOLUTION? 

3.1 WZ Sge as a magnetic propeller 

In sections l2.ll and f2.2l we have argued that the value of a co id 
in WZ Sge may not be different from that in other DN. This 
is consistent with the observed AM aC c, if -Rout ~ 0.5a, and 
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Figure 8. Surface density profiles of the disc in WZ Sge as simulated in the one-dimensional code with a magnetic torque applied. Four 
snapshots are shown in frames (a) to (d). Frame (a) shows the disc just before the onset of an outburst. Frame (b) is seven hours later 
showing the heating waves propagating inwards and outwards. Frame (c) shows the disc four days later, at the height of the outburst, 
with the disc in the fully hot state and with no magnetic truncation. Finally, frame (d) shows the early inter-outburst period, eighteen 
days later, with a greatly depleted disc. The critical trigger densities are plotted with dashed lines. 



with the decay time of outbursts. However the recurrence 
time problem (Smak 1) still remains unresolved. 

The observation of a 27 .87 oscillation in the quiescent 
x-ray lightcurve of WZ Sge (Patterson 1998) suggests that 
the white dwarf is magnetic. The introduction of a magnetic 
field to a rapidly rotating white dwarf raises the possibility of 
increasing the recurrence time by removing the inner, most 
unstable regions of the accretion disc. This is a result of the 



magnetic torque which can, rather than accreting the inner 
disc, force mass close to the white dwarf to larger radii. This 
magnetic propeller effect occurs when the stellar magnetic 
field moves more rapidly than the disc material so that the 
disc gains angular momentum at the expense of the star. 

The inner disc can be depleted by the propeller and the 
surface density profile of the disc will then be altered as ma- 
terial is forced to orbit further out in the disc, at radii close 
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Figure 9. Grey-scale density plots of the accretion disc in WZ Sge, simulated in the two-dimensional SPH code. The first plot shows 
the disc in the non-magnetic case after 79 orbits. The second and third plots show the disc in the magnetic case, after 45 and 79 orbits 
respectively. The first of these two plots is just before the disc begins to become eccentric and the second shows the later, eccentric disc. 
The primary Roche lobe is plotted as a solid line and the primary star is marked by a circle (not to scale) in all three plots. 
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Figure 10. Azimuthally averaged plots of surface density against radius in two dimensional simulations of WZ Sge. The first plot shows 
the disc profile in the magnetic and non-magnetic cases after 45 orbits, before the disc becomes eccentric. The second plot shows the 
disc after 79 orbits, when the disc has become eccentric in the magnetic case. Note that the £ axis scale is not the same in both plots. 



to the tidal limit. The magnetic propeller inhibits accretion 
so that mass is forced to build up in the outer disc, where 
the critical surface density is higher and the E evolution 
time-scale (E cr it/S) is longer. The magnetic propeller effect 
cannot therefore only increase t IOC but can also allow the disc 
to accumulate enough mass prior to outburst to supply an 
increased AM aC c, with the outer disc acting as a large reser- 
voir of mass. In this way both (Smak 1) and (Smak 2) 
can be explained simultaneously with a standard value for 
a co id, and without appealing to episodes of enhanced mass 
transfer. 

It should be noted that disc truncation is not normally, 
in itse lf, sufficient to increase recurrence times. iMenou et al.l 
(2000) show this by considering the evaporation of the inner 
disc. The fundamental difference between the magnetic pro- 
peller and other mechanisms such as evaporation and mag- 
netically enhanced accretion is that the propeller causes a 
build-up of mass in the outer disc by preventing accretion 
onto the white dwarf (during quiescence only). The overall 
profile of the disc is therefore completely altered, permitting 
much more mass to be stored in the disc. The mechanism 
presented here is analogous to that which was proposed by 



iMatthews et alJ (|2004j ) to explain long recurrence times in 
FU Ori type young stellar objects. 

The magnetic propeller should, in forcing the disc out- 
wards, aid access to the 3:1 resonance. It is likely there- 
fore that fur ther mass wi ll be stored in the outer disc, as 
proposed by lOsakil (|l995l ). providing an even larger reser- 
voir AM for mass between outbursts. This tidal effect may 
therefore amplify the result of the weak propeller driving 
recurrence times to become even longer. 

In the remainder of this section we outline the theoret- 
ical basis for the magnetic propeller model, and present the 
results of numerical experiments. 



3.2 The magnetic model 

We adopt a description for the magnetic interaction which 
assumes that as material moves through the magnetosphere 
it interacts with the local magnetic field via a velocity de- 
pendent acceleration of the general form: 

dmag = —k (V - Vf) ± (8) 
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Figure 11. Evolution of the peak surface density in the disc 
of WZ Sge as a fraction of the critical surface density in non- 
magnetic and magnetic cases as a function of time. Results are 
taken from two-dimensional SPH calculations. 



Figure 6. Plot showing accretion onto the white dwarf as a func- 
tion of time in magnetic and non-magnetic cases on the same time 
scale as Fig [5] The magnetic outbursts have a higher peak accre- 
tion rate than the standard outbursts. 
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Figure 7. Plot showing viscous dissipation for the case of a mag- 
netically truncated disc close to marginal stability with B ~ 1 kG. 
The recurrence time is increased to t rcc ~ 20 yr. 



where v and v/ are the velocities of the material and mag- 
netic field respectively, and the suffix _L refers to the veloc- 
ity components perpendicular to the field lines. A version of 
this description has been succe ssfully appl ied to a number 
of other intermediate polars bv lKinel (ll99St): IWvnn fc Kind 
l|l995h : IWvnn et all (|l997r i and iKing fc Wvnnl (|l999T i. It 
has also been a pplied to d iscs around you n g sta rs by 
iMatthews et al.l 1|2004 ) and iMatthews et al.l i|2005l ) The 



magnetic acceleration in equation (|8]) is intended to repre- 
sent the dominant term of the magnetic interaction, with k 
playing the role of a 'magnetic a'. The magnetic time-scale 
can be written in terms of k as 



ft" 1 , 



Vf 



(9) 



It is assumed that the magnetic field of the WD is dipolar 
and that it rotates as a solid body with the star. We refer to 
a system as a weak magnetic propeller when, in quiescence, 
the magnetic field is strong enough for the system to behave 
as a magnetic propeller, but not strong enough to eject a sig- 
nificant fraction of the transferred mass from the Roche lobe 
of the primary (as is the case in AE Aqr). During outburst 
the field should not be strong enough to prevent accretion. 
This requires the following hierarchy of time-scales: 



(10) 



where £ v i sc ,h 

scales respectively. Assuming the temperature of the quies- 
cent disc to be T ~ 5000 K and a c = 0.01, and that in 
outburst we have T ~ 5 x 10 4 K and oth~ 0.1, equation {(Sj 
gives 



1 days < fmag < 40 days. 



(11) 



A further condition for a system to act as a weak propeller 
is that the magnetic interaction takes place at radii exceed- 
ing the co-rotation radius _R co = (GMiP S p in /47r 2 ) 1,/2 , where 
the fieldlines rotate with the local (circular) Kepler veloc- 
ity. A spin period of 27.87 s for the WD in WZ Sge gives 
Rco > Rwd ~ 10 9 cm, easily fulfilling this requirement. 
If, conversely, the interaction were to take place within R co 
then the effect of the magnetic field would be to enhance 
accretion. 

It now remains to quantify i mag . There are a number 
of different models for the interaction of the inner disc and 
the white dwarf magnetosphere which have been applied to 
magnetic CVs. Here we examine the two extreme cases in 
which the accretion disc is completely magnetized, and con- 
versely, in which the disc is diamagnetic. In both cases we 
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still assume that the stellar field has a dipolar geometry, and 
that any local field distortions may be treated as perturba- 
tions on this structure. 

When the magnetic field lines permeate the inner disc 
any vertical velocity s hear tends to amplify the toroidal field 
component B v (e.g. lYi fc Kenvonl [l997l ) . Diffusive losses 
counteract this amplification and in steady state, these ef- 
fects balance to give 



(12) 



Here 7 is a parameter that accounts for the uncertainty in 
the vertical velocity shear and is of the order unity, r < 1 rep- 
resents an uncertainty in the diffusive loss time-scale, whilst 
and J1k are the angular velocities of the WD and disc re- 
spectively. The tension in the field lines results in a magnetic 
acceleration which may be approximated by 



B v 




7(n*-n K ) 









1 

PdTc 



4% 



(13) 



where pd is the disc density, n is the unit vector perpendic- 
ular to the field lines and r c is the radius of curvature of 
the field lines. It is possible to parametrize r c as 



(14) 



where e is a constant of order unity. Combining equations 
(fl"2|) . (fT4|) and (fl"3|) . with H ~ 0.17?, the magnetic accelera- 
tion is given by 



5~/B z 



2-irRrepdVK 
and i mag by 

2TvRrep d VK 



(v K - v f ) 



■!'K 



5jB z 



\VK — V{\ ± 



(15) 



(16) 



At the circularisation radius i? c irc ~ 10 10 cm, the disc den- 
sity is pd ~ 10 -6 g cm -3 (from equation (J3]l, with H ~ 
0.1-R). Assuming r ~ t ~ 7 ~ 1 we are able to constrain the 
magnetic moment p ~ B Z R 3 using (|ll|l . giving 



2 x 10 32 G cm 3 < n < 10 33 G cm 3 



(17) 



For the opposite case of a diamagnetic flow we assume 
that the inner regions of the disc are broken up into a se- 
ries of individua l filam ents, as predicated by the analysis of 
lAlv fc Kuijpersl (|l99Ch . In order to follow the evolution of 
the diamagnetic gas blobs (which are formed we ll outside the 
co-rot ation radius) we follow the approach of IWvnn et al.l 
1 19971 ). This approach assumes that the diamagnetic blobs 
interact with the field via a surface drag force, characterized 
by the time-scale: 



CAPblbB 



r<'K 



(18) 



\vk - Vf\± 

where ca is the Alfven speed in the medium surrounding 
the plasma, B is the local magnetic field, pb is the plasma 
density, and lb is the typical length scale over which field 
lines are distorted. Making the approximations that pth ~ 
£ C rit(10 10 cm), ca — c s and v ~ vk we find that 

5 x 10 30 G cm 3 < n < 3 x 10 31 G cm 3 . (19) 



3.3 One-dimensional numerical results 

In this section we use a one-dimensional numerical model 
to confirm that it is feasible to create a depleted region in 
the centre of a quiescent accretion disc with a propeller. 
We also show that the propeller can lead to a build up of 
mass in the outer disc, and to extended recurrence times. 
The code o utlined below is based on that used to simi- 
lar effect in iMatthews et all |2004) . but with the addition 
of viscosity switching to reproduce the effect of disc out- 
bursts. A thin axisymmetric accretion disc is assumed to 
be rotating about a central mass M. The hydrodynamical 
equations are averaged over the azimuthal dimension so that 
non-axisymmetric effects, such as tidal forces from the sec- 
ondary star, are not modelled. It is assumed that the disc is 
sufficiently cold that the dynamical time-scale t^yn ~ R/v v 
is much smaller than the viscous or magnetic timescales so 
that a Keplerian approximation can then be adopted for the 
azimuthal motion. The magnetic torque due to the rotating 
WD is added using a parametrization of the form 

1 = ^f* 1/2 , (20) 



A = — 

t A 

where G represents the universal constant of gravitation, 
M is the mass of the central star and t/i is the time-scale 
on which the local disc material gains angular momentum, 
which in this case is given by equation (|16p . The resulting 
equatio n for the evolution of s urface density with time (de- 
rived in IMatthews et al.ll20o4 ) is 



dE _ 3__d_ 

dt RdR 

_ (3__d_ 
RdR 



R 



1/2 



d{R l/2 uE) 



1 



R 7/2 

where j3 is defined as 



OR 
R 

Rco 



(21) 



3/2 



A* 



(22) 



2ttVGM 

We adopt the lShakura &: Sunvaevl (|l973l ) viscosity prescrip- 
tion © and the viscous dissipation is calculated using equa- 
tion ©. The thermal viscous disc instability (DIM) must 
be modelled if disc outburst behaviour is to be reproduced, 
and the viscosity is therefore permitted also to vary in 
time. The approach to viscosi ty switching adopt ed here is 
very similar to that used by (Truss et al.l (|2000h in their 
smoothed particle hydrodynamic simulations. Disc annuli 
are switched between high and low viscosity states, by al- 
tering the_viscous_0£parMieter, using critical density 'trig- 
gers' (C annizzo et al.lll988l ) which lie at the extremes of the 
limit cycle. Accordingly when the local surface density E 
exceeds that of equation the local viscous parameter is 
altered smoothly on a thermal timescale (tth ~ ( a ^) _1 ) to 
a value of Qhot = 10a co id = 0.1. Later, when E falls below 
a lower trigger density, the viscous parameter is returned to 
its original value a co id = 0.01. In the high state the accretion 
disc should also be hotter, so the temperature T is switched 
between 5000 K and 50000 K in the same way. 

The code calculates the diffusive and advective terms 
separately on a grid of resolution 500 and combines the re- 
sults using operator splitting. The diffusive term is calcu- 
lated using a forward in time space-centred scheme (FTSC) 
and the values of S produced here are used as input for a 
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modified two-step Lax-Wendroff scheme l|Press et al.lll992h 
which is used to solve the advective term. The inner bound- 
ary, representing the stellar boundary layer, is placed at 
R+ = 1.0 x 10 9 cm and an outflow condition such that ^-M- 
is constant is applied. An inflow condition is applied at the 
outer boundary which is placed at the 0.4a, following SR. 
The stellar mass is set to M = 1.2 Mq and mass is added 
at the outer boundary at a rate of M — 2 x 10 -10 Mq yr -1 . 
The white dwarf spin period is set to P sp i n = 28 s which 
gives a co-rotation radius of R co = 1.5 x 10 9 cm. 

Calculations were performed without a primary mag- 
netic field and with a field of B z ~ 1 kG (p ~ 10 30 Gem 3 ), 
although the precise value depends on the prescription 
adopted for the magnetic acceleration. The effect of the mag- 
netic field was to truncate the inner disc, producing a de- 
pleted central region and to reduce the peak surface density. 
The inner disc is likely to be the starting point of any out- 
burst in the non-magnetic case while in the magnetically 
truncated case this should instead occur close to the trun- 
cation radius. It is important that in neither case is the disc 
close to marginal stability, where the steady state profile 
would peak very close to E/E cr i t = 1. Such marginal stabil- 
ity is another possible reason for long outburst times, but it 
would not cause an increase in outburst amplitude. 

When the viscosity switching mechanism was enabled 
a series of outbursts occurred in both the magnetic and 
non- magnetic cases. These are plotted in Fig [3] In the non- 
magnetic case the outburst time is i Q b ~ 10 d and the re- 
currence time is t rcc ~ 20 d. This is roughly in agreement 
with observations of normal DN, although the outbursts are 
longer than expected. The outbursts vary in amplitude but 
are of the expected order of magnitude. In the magnetic case 
the recurrence time is increased by ~ 40 times to t rcc ~ 2 yr, 
and the outburst amplitude is increased by a factor of ~ 2. In 
both cases the quiescent dissipation is more than two orders 
of magnitude less than the outburst dissipation. The profile 
of an individual magnetically moderated outburst is shown 
in Fig [4] The outburst lasts for toh ~ 10 d, notably shorter 
than is observed. It appears, after an initial sharp rise, to 
rise on a time-scale of a few days. This is close to the hot vis- 
cous timescale, rather than to the thermal time-scale. This 
is probably because once the disc becomes hot, on a thermal 
time-scale, mass migrates into the centre of the disc, filling 
the evacuated central regions on the viscous timescale. The 
outburst decays on a similar viscous time-scale, and then 
much more quickly, on the thermal time-scale, as the disc 
returns to the cool state. 

Fig [5] shows how the mass of the accretion disc varied 
during outbursts. The total peak disc mass in the magnetic 
case was Mdi sc ~ 8 x 10 23 g, twice that of the non-magnetic 
case. Further, the mass consumed in the magnetic outbursts 
was A ~ 4 x 10 23 g, ten times that consumed in the non- 
magnetic outbursts, half of the total accretion disc mass, 
and approaching the value deduced from observations. We 
also note that the mass growth rate in the magnetic case 
is close to linear, this confirms that the disc is not close to 
marginal stability, as was discussed above. The results from 
the non-magnetic case show a series of alternate outbursts 
and mini-outbursts. Mass accretion onto the white dwarf 
is plotted in Fig [SJ on the same time axis as Fig [S] The 
accretion rate peaks at 5 x 10 -9 Mq yr -1 in the magnetic 
case, which is a factor of 2 more than the non-magnetic case. 



It is also notable that the cold viscous decay, which is not 
visible in Fig [4] can be clearly seen in the magnetic case. 

In Fig [5] the detailed evolution of a magnetically mod- 
erated outburst can be followed. During quiescence the pro- 
peller has prevented the disc from spreading in to the star, 
and by preventing accretion has forced a large build up of 
mass in the outer disc. The outburst begins when the surface 
density first exceeds the critical value, close to the trun- 
cation radius, rather than near the boundary layer in the 
usual way. Heating waves then sweep inwards and outwards 
on a thermal time-scale until the whole disc is in the hot 
state. The viscosity of the disc is now higher, and the vis- 
cous timescale shorter, allowing mass to spread inwards and 
quench the propeller. The disc now takes a form close to the 
standard non-truncated Shakura-Sunyaev solution. Mass ac- 
cretion onto the white dwarf is high and the density of the 
disc drops rapidly until the lower trigger is reached at the 
outer edge of the disc. At this point a cooling wave sweeps in- 
wards so that the disc becomes cool again and is gradually 
shifted outwards as the propeller is reestablished. Around 
half of the total disc mass is consumed in the outburst. 

Results from these calculations confirm that a magnetic 
propeller is capable of truncating the inner accretion disc 
and may also lead to an extended recurrence time. However 
the recurrence time is still an order of magnitude shorter 
than that observed with these simulation parameters. Two 
possible mechanisms for further extending the recurrence 
time suggest themselves. The first requires only a slightly 
stronger magnetic field. If the magnetic field is stronger 
then the disc will be truncated at a larger radius, and this 
can be tuned so that the disc is very close to the steady 
state when it reaches the outburst trigger. As the disc ap- 
proaches this steady state its growth rate slows asymptoti- 
cally to zero so that the recurrence time in th is case can be 
extremely long. This effect was discussed by IWarner et al.l 
l|l996t h iMenou et all (|200Cj ). It is distinct from, and acts 
in addition to the mechanism described above which is con- 
cerned instead with the linear regime of disc growth. Results 
from a calculation using this marginal stability effect are 
shown Fig. [7] The recurrence time is increased to t lac ~ 20 yr 
but neither the mass consumed nor the amplitude of the out- 
bursts is greatly altered from the standard magnetic case il- 
lustrated in Fig|3j so that the marginal stability effect would 
be unable, in itself, to solve (Smak 2). This mechanism 
also requires fine tuning; if the field is slightly too strong 
outbursts will be prevented entirely. 

A second possible explanation for the further extension 
of the recurrence time, which does not require such fine 
tuning, is that the size of the mass reservoir may be fur- 
ther increased by tidal effects. The one-dimensional treat- 
ment detailed above precludes the consideration of binary 
tidal forces and resonances. To better model the structure 
of the accretion disc, including its non-axisymmetric modes, 
we turn to the complementary two-dimensional model de- 
scribed in section 13.41 

3.4 Two-dimensional numerical results 

Further numerical calculations were performed in two 
dimensions with smo othed particle hydrodynamics 
(SPH) simulations (e.g. iKunze. Speith fc Hessmanl 1200 ll ; 
ISchafer. Speith fc Hippl 120041 ). The above prescription {T5} 
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for the magnetic acceleration was implemented in an SPH 
code and calculations were performed with and without a 
dipolar magnetic field in the plane of the accretion disc. 
The system parameters were similar those used in the one- 
dimensional simulations described above, although the mass 
transfer rate was M — 3 x 1CP 11 M© yr -1 . Each calculation 
used around 40,000 particles with a fixed smoothing length. 
Particles were injected at the first Lagrangian point and 
removed when they approached the primary star. 

The resulting discs are shown in Figs [9] and 1101 plot- 
ted as face-on disc views and radial profiles respectively. 
The profiles are somewhat different to those from the one- 
dimensional simulations, shown in Fig [5] as a result of the 
torques from the binary potential. The effect of the magnetic 
field in pushing the disc outwards is more pronounced than 
in the one dimensional model. In the magnetic case, because 
matter is driven towards the resonant radii, the disc becomes 
eccentric and begins to develop some spiral structure. These 
are clearly seen as two peaks in the second plot of Fig 1101 

This structure also affects the time evolution of the sur- 
face density. Fig 1111 shows how the peak surface density as 
a fraction of the critical surface density grows with time in 
the magnetic and non-magnetic cases. This can be viewed 
as a measure of how close the disc is to outburst, since when 
any substantial part of the disc reaches its local critical den- 
sity an outburst will be triggered. In the non-magnetic case 
there is a simple monotonic growth of this parameter. The 
noise visible in the otherwise very smooth trend is caused 
by density enhancements associated with accretion events 
at the inner boundary, which may be a result of the limited 
resolution of the simulation. This monotonic growth is as 
expected from Fig [10] which shows that the non-magnetic 
disc profile keeps the same form, with the density peak re- 
maining at the same point around ~ 0.3a. If surface density 
continues to grow at this rate until outburst then this sim- 
ulation represents around 3% of the duty cycle for a non- 
magnetic system. In the magnetic case the same parameter 
grows monotonically at first, much more slowly than in the 
non-magnetic case. Its behaviour becomes stochastic once 
the disc is driven to become eccentric. It is problematic to 
extrapolate from such a short period to the length of the 
recurrence time, however the overall trend remains clear: 
the approach to outburst is clearly slower in the magnetic 
case. This result does not contradict that of 13.31 It therefore 
seems unlikely that effects due to the binary would prevent 
the propeller from leading to a long recurrence time. It also 
remains possible that the tidal forces could lead to the de- 
velopment of a much larger reservoir, and hence a longer 
recurrence time, than in the simple axisymmetric case. 



4 OBSERVATIONAL CONSEQUENCES 

The weak magnetic propeller model for the outbursts of 
WZ Sge results in a number of observationally testable pre- 
dictions. The major consequence of the model, the exis- 
tence of a substan tial hole in the inner d isc, h as already 
been suggested by iMennickent fc Arenas! ([l998) . The au- 
thors estimate the ratio of the inner and outer disc radii 
(R — Rin/Rout) from observations of the H a emission line 
widths. WZ Sge was found to show an extremely large value 
of 7? ~ 0.3. This value agrees reasonably well with the sec- 



ond plot in Fig [9] where Rout — 0.5a and Ri n ~ 0.2a for the 
magnetic case, and is much larger than the value expected 
if the accretion disc extended all the way to the WD surface 

(R~om). 

The asynchronous magnetic cataclysmic variables (the 
intermediate polars) show multiple periodicities in their light 
curves. WZ Sge has been observed to display oscillations at 
27.87 s and 2 8.96 s as well as we aker oscillations at 28.2 s 
and 29.69 s ( Lasota et al.l Il999l . and references therein). 
lLasota et al.l 1 19991 ) interpret WZ Sge as an intermediate 
polar, and identify the 27.87 s oscillation as arising from 
the spinning white dwarf, and the other oscillations as beat 
periods between the white dwarf spin and material orbiting 
within the disc. The beat periods arise from the reprocess- 
ing of x-ray emission from the WD on the surface of the 
disc, the hotspot and/or the secondary star. Assuming the 
material in the disc is moving in circular, Keplerian orbits 
the radius associated with these beat periods Rbcat can be 
written as 

p£gm\ 1/3 

4tt 2 J 



R 



(23) 



where 
1 1 



K 



1 

-fobs 



(24) 



Here Pk is the Keplerian period of the disc material, and 
-Fobs is the observed beat period. Assuming M — 1.2Mq we 
find i?28.96 ^ 1.3 x 10 10 cm (~ 0.3a), R 28 .2 ^ 2.8 x 10 10 
cm (~ 0.6a), and i? 2 9.69 9.4 x 10 9 cm (~ 0.2a). The radii 
F?28.96 and -R29.69 correlate reasonably well with material 
orbiting close to the inner edge of the quiescent disc in two- 
dimensional numerical models (cf Fig O - The radius R2S.2 
on the other hand, is quite close to the hotspot (stream/disc 
impact region) which should be located at ~ 0.5a, as seen in 
the magnetic cases in Fig [5] Obtaining a closer fit between 
these radii and the model would require tighter constraints 
on the WD mass, and q, as well a suitable model for the x- 
ray reprocessing within the disc. Nonetheless absence of beat 
radii from the supposed evacuated region is encouraging for 
the magnetic propeller model. The fact that the oscillations 
are not observed during outburst is also entirely consistent 
with this model: the magnetosphere is compressed during 
outburst, and the system is observationally "non-magnetic". 
The model outlined in this paper predicts that WZ Sge has 
an eccentric precessing disc, and it might therefore be ex- 
pected to display permanent superhumps. 



5 THE SPIN EVOLUTION OF THE WHITE 
DWARF 

The spin evolution of the white dwarf in WZ Sge is de- 
termined by angular momentum transfer in the quiescent 
and outburst phases. In quiescence the WD acts as a mag- 
netic propeller, i.e. material approaching the WD is forced 
to larger radii. We can approximate this interaction by as- 
suming the transferred gas, with initial specific angular mo- 
mentum (GM Rdrc) 1 ^ 2 , is forced to orbit at a mean outer 
radius -Rout- The associated spin down torque on the WD is 
then 



J d = -CM 2 (GM) 1/2 [Rl^-R. 



,1/2 



(25) 
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where the mass transfer rate accross the LI point Mi ~ 
10 15 gs _1 , .R ou t ~ 0.5a, i? c irc ~ 0.35a, and the parameter 
( > 1 if material at flout loses angular momentum to the 
secondary star via tides (this will be expected to occur once 
the outer edge of the disc reaches the tidal radius). Con- 
versely, in outburst, the WD gains angular momentum with 
accreted matter at the rate 

J u = M acc (GMR in ) 1/2 , (26) 

where we ad opt an i nner disc radius of Ri n = 10 9 cm. From 
the results of lSmakl <| 19931 ) the average accretion rate during 
outburst can be estimated as M acc ~ 10 18 g s _1 . 
The evolution time-scale of the white dwarf spin period can 
be written in terms of J as 

P spin _ 2tvI 



p . 

1 spin 



P ■ J 
- 1 spiiW 



(27) 



where I is the moment of inertia of the WD. Averaging over 
a duty cycle the torque on the WD is given by 



< J >: 



(28) 



(ircc ~\~ £o) 

where t b ~ 50 d is the outburst duration and we take t rcc ~ 
30 yr. From the estimates above (assuming £ ~ 1) we find 
Juto/ Jdtrcc ~ 6. The greater spin-up torque gives a spin up 
on the time-scale 



2-kI 



Td 



| P\ spin | -fspin I < J > 



10 9 yr. 



(29) 



The WD in a weak propeller may experience a net spin up, 
as above, or a net spin down, if outbursts are shorter or the 
mass transfer rate is higher for example. Both these results 
may lead to a spin up/spin down cycle. 

Net spin up: In this case, illustrated by WZ Sge, the 
net spin evolution is towards shorter periods. Therefore, 
according to equation © the propeller will become stronger 
with time. The condition t mag < t v isc,h will eventually 
be satisfied. At this point accretion will be prevented 
completely and outbursts will cease; the system will become 
AE Aqr like. This phase should be short lived however and 
it must lead to a spin down. Therefore a cycle between 
WZ Sge and AE Aqr type stars may arise. 

Net spin down: In this opposite case, the net evolution 
is towards longer white dwarf spin periods. The propeller 
becomes weaker with time until it is eventually quenched 



when i n 



^ tvis 



The system would then look like a nor- 



mal SU UMa star. This SU UMa phase would usually result 
in white dwarf spin up. Therefore an alternative cycle may 
exist between WZ Sge and SU UMa type systems. Both of 
these cycles are a natural consequence of hierarchy (|10[) . 



6 DISCUSSION 

We have argued that the outbursts of WZ Sge can be ex- 
plained using the standard DIM value for a co id if we as- 
sume that the white dwarf in the system acts as a magnetic 
propeller. Numerical models of the magnetic disc, in both 
one and two dimensions, predict an increased inter-outburst 
time free, and a disc massive enough to fuel an increased 



outburst accretion rate. No episode of enhanced mass trans- 
fer from the secondary star is required to trigger outbursts, 
or to supply mass during outburst. 

The weak-propeller models assume that the magnetic 
tension force is the dominant cause of angular momentum 
transfer in the quiescent disc, and that the force is propor- 
tional to the local shear between the disc plasma and the 
magnetic field. In outburst angular momentum transfer is 
dominated by viscous diffusion, and mass accretion on to 
the white dwarf is permitted. This results in an estimated 
WD magnetic moment in the range 



2 x 10 32 G cm 3 < n < 10 33 G cm 3 



(30) 



if the disc is assumed to be fully magnetized, or a value in 
the range 



5 x 10 30 G cm 3 < fj, < 3 x 10 31 G cm 3 



(31) 



if the disc is assumed to be diamagnetic. These values for the 
WD magnetic moment would suggest that WZ Sge is a short 
period equivalent of the intermediate polars (IPs). Most IPs 
lie above period gap (P or b > 3 hours) and have magnetic 
moments in the ran ge 10 31 G cm 3 < (j, < 10 34 G cm 3 



iKing fc Wvrml (|l999h pointed out the only two confirmed 
IPs with Porb < 2 hours (EX Hya and RX1238-38) have 
A* ?1 10 33 G cm 3 and long spin periods (67 min and 36 min 
respectively). WZ Sge would therefore be the first weakly 
magnetic /i < 10 32 G cm 3 CV to be found below the period 
gap. 

In differentiating between the models of the disc- 
magnetosphere interaction it is useful to consider the case of 
AE Aquarii. This system is thought to contain a WD with a 
magnetic moment ~ 10 32 G cm 3 , which ejects ~ 99% of th e 
mass transfer stream from the system l|Wvnn et al.lll997h . 
If the WD magnetic moment (j, > 10 31 in WZ Sge then it 
would satisfy the criteria necessary to become an ejector 
system. Hence in the case a fully magnetized disc, it would 
seem that the presence of the disc itself is the only protection 
the system has from ejecting the transferred mass. If the disc 
was completely accreted or destroyed for any reason WZ Sge 
would resemble a short period version of AE Aqr. This would 
not apply to the case where the WD has a magnetic moment 
in the range predicted by the diamagnetic model. 

A WD spin evolution cycle results as a natural con- 
sequence of the time-scale hierarchy (|10p . In the case of 
WZ Sge we expect the spin cycle to run between the current 
state, characterised by long inter-outburst times, and that 
of an AE Aqr-like strong propeller. In some similar systems 
however, spin-up phases would look like a normal SU UMa 
system, whilst the spin down phases would be WZ Sge like. 
WZ Sge-like systems following such spin cycles are expected 
to be relatively rare because of the constraint (|10p . which 
results in the restricted range of allowed WD magnetic mo- 
me nts presented above. 

iMennickent fc Arenas! (| 1998f ) find evidence for ring-like 
accretion discs in long supercycle length SU UMa stars from 
a radial velocity study of the H a emission lines of these sys- 
tems. WZ Sge is the most extreme of these objects with a re- 
currence time of ~ 12, 000 d and R = R^/ Rout = 0.3, which 
is in good agreement with our numerical results. Other 
SU UMas show a strong positive correlation between su- 
percycle times and the R parameter. This result is in agree- 
ment with the model presented in this paper for WZ Sge. 
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The WDs in the other, less extreme, SU UMa stars would 
simply act as less efficient propellers: the WDs would have 
weaker magnetic moments, or lower spin rates (which may 
reflect various stages of the spin up/down cycles postu- 
lated above). Strong candidates for WZ Sge-like systems 
are RZ Leo (t rcc < 4259 days, R ~ 0.16), CU Vel (i rcc ~ 
700 - 900 days, R ~ 0.1 5) and WX Cet (t rec ^ 1000 days, 
R ~ 0.12), (taken from iMennickent fc Arenas! Il998l , table 
5). 
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